Intersite coupling effects in a Kondo lattice 
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The La dilution of the Kondo lattice CeCoIns is studied. The scaling laws found for the magnetic 
susceptibility and the specific heat reveal two well-separated energy scales, corresponding to the 
single impurity Kondo temperature Tk and an intersite spin- liquid temperature T* . The Ce-dilute 
alloy has the expected Fermi liquid ground state, while the specific heat and resistivity in the 
dense Kondo regime exhibit non-Fermi-liquid behavior, which scales with T* . These observations 
indicate that the screening of the magnetic moments in the lattice involves antiferromagnetic intersite 
correlations with a larger energy scale in comparison with the Kondo impurity case. 



PACS numbers: 71.27.+a, 75.20.Hr, 75.30.Mb 

The nonuniversal behavior of heavy fermion com- 
pounds has challenged experimental and theoretical 
physicists over many years. A wide range of low- 
temperature phenomena, such as strongly enhanced 
paramagnetism, magnetic/quadrupolar long-range order, 
Kondo insulators, unconventional superconductivity, and 
non-Fermi-liquid behavior, have been reported. This va- 
riety of phenomena is believed to arise through the com- 
petition and interplay of the Kondo effect, band struc- 
tures and intersite correlations. Both, the Kondo-esque 
resonance at the Fermi level JjJ and the Ruderman-Kittel- 
Kasuya-Yosida (RKKY) interaction in a Kondo lattice, 
are the consequence of the large Coulomb repulsion in the 
partially filled 4/(5/)-shell in conjunction with the hy- 
bridization of the /-electrons with the conduction state. 

In the single impurity limit, the low energy transport 
and the thermodynamic properties scale with a single 
parameter, the Kondo temperature Tk B. As T is de- 
creased below Tk, the moments which are localized at 
high temperatures are screened into local singlets by the 
conduction electrons, giving rise to a local Fermi liquid 
state. The agreement between experiment and theory is 
excellent for many dilute systems 0] . 

For Kondo lattice systems, it has been a long-standing 
issue how the intersite interactions affect the local char- 
acter of the Kondo screening. If the singlet formation in 
the lattice involves the neighboring localized moments, 
it is expected that the spin screening is accompanied by 
antiferromagnetic (AF) short-range correlations, which 
may modify Tk or yield a new energy scale. While a 
study of the Kondo lattice within the Gutzwillcr approx- 
imation || predicts such a "lattice enhancement" of Tk , 
1/N expansions Q] indicate local singlet formation with 
no change in the energy scale. The lattice enhancement 
has been confirmed by Monte Carlo, exact diagonaliza- 
tion, and density matrix renormalization-group H for the 
1 D Kondo lattice, but backscattering across the Fermi 
surface may affect this result. Monte Carlo studies for 
the 2D half-filled Kondo lattice yield different renormal- 
izations of the spin and charge gaps B . 



Experimentally, the strength of the intersite coupling 
can be varied by substituting the magnetic Ce ions by 
their nonmagnetic analog La. This approach continu- 
ously interpolates between the Kondo lattice and the 
single impurity limit. Several systems have been stud- 
ied, e.g. (Ce,La)Al 2 , (Ce,La)B 6 , and (Ce,La)Cu 6 0, but 
it is hard to arrive at quantitative conclusions, since La 
substitution often changes both the crystalline electric 
field (CEF) parameters and Kondo coupling. 

The ideal conditions of concentration independent Tk 
and CEF are met for (Ce,La)Pb3 ||. It was reported 
that the specific heat and susceptibility scale with the Ce- 
concentration for the entire alloying range and follow the 
predictions of the 5=1/2 Kondo impurity model with 
Tk — 3 K ||. No signs of intersite couplings were found, 
except in the concentrated limit, where they eventually 
lead to the AF ordering at T N ~ 1 K for CePb 3 . 

In this Letter we present experimental evidence for 
three well-separated energy scales in Cei-zLazCoIr^: 
Tk, the energy scale for intersite interactions T*, and the 
level splitting due to CEF. CeCoIns is a heavy fermion 
unconventional superconductor with a high transition 
temperature T sc = 2.3 K With La dilution, 

Cei_a;La K CoIn5 has a paramagnetic ground state with 
Ce keeping essentially the same valence and CEF split- 
ting through the entire alloying range. We found scaling 
functions with the energy scale T* for the deviations of 
the specific heat and the susceptibility from the single im- 
purity behavior. This scale T* is attributed to intersite 
interactions. The specific heat in the dense Kondo regime 
exhibits non- Fermi-liquid behavior with the characteris- 
tic temperature T*, while for dilute alloys it follows that 
of the S — 1/2 Kondo impurity model. The deviation 
of the susceptibility from the single ion behavior shows 
a strong decrease below T* , suggesting the formation of 
a spin-liquid. The scaling procedure reveals that T* sys- 
tematically increases from the single-ion Tk with increas- 
ing Ce concentration. This indicates that the screening 
process of the localized moments in the lattice has higher 
energy scale, involving AF intersite correlations. 
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FIG. 1. Magnetic in-plane resistivity p m for various x of 
Cei-^La^CoIns. Right inset: The log-log plot for the inelastic 
part of p m vs. r/T co h; the solid line is the T linear fit and the 
vertical lines mark the onset of superconductivity. Left inset: 
pm vs. T 2 for the incoherent regime. 

Single crystals were prepared using an In self-flux 
method {J[l2). Before growing the crystals, Ce and La 
metals are premixed by arcmelting at least 10 times to get 
a homogeneous mixture. X-ray powder diffraction pat- 
terns show a single phase tetragonal HoCoGas structure, 
which can be viewed as alternating layers of (Ce,La)In3 
with layers of CoIn2 along the c axis. The resistivity was 
measured by standard four-probe dc and ac techniques 
down to 25 mK using the low temperature facilities at the 
NHMFL. The magnetization was measured down to 1.8 
K with a Quantum Design MPMS SQUID magnetome- 
ter and the specific heat (Cp) by a thermal relaxation 
method using a Quantum Design PPMS. 

Figure 1 shows the magnetic part of the in-plane 
resistivity p m defined as the difference in p(T) of 
(Ce,La)CoIn5 and LaCoIn5 (the nonmagnetic, isostruc- 
tural, analog of CeCoIn5) divided by the Ce concentra- 
tion, 1— x. Above 50 K, the p m data for all x collapse onto 
a single curve with the —log T dependence characteristic 
of the Kondo effect. Hence, the high T Kondo tempera- 
ture, T£, is essentially independent of a;. In the Ce rich 
region, p m shows coherent metallic behavior after rising 
through a peak that defines the coherence temperature 
Tcoh- The La dilution reduces T co h, which tends to zero 
near x = 0.5, where p m saturates below 500 mK. With 
further La dilution, p m evolves toward single ion behavior 
with a log T dependence down to 500 mK that suggests 
a small Tk of about 1-2 K. Hence, x — 0.5 roughly sepa- 
rates the coherent and single ion regimes. It is interesting 
to note that 50 % is close to the percolation limit of 41 
% of a 2D square lattice (consistent with the 2D nature 
of Ce-Ce nearest neighbor network in the Celn3 layers). 
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FIG. 2. Excess susceptibility A\(T) = X (T) - Xo{T) 
(xo(T) refers to the single impurity limit at x = 0.955), scaled 
with T/T* . The solid line corresponds to an activation type 
fit. Inset: the c-axis component of M/H versus logT. 

The inset of Fig. 2 shows the c-axis component of the 
susceptibility Xc{T) per mole Ce vs. logT. As for p m , 
the data above 80 K collapse onto the same curve for all 
concentrations, again indicating that T^i is independent 
of x. The high T data follows the modified Curie- Weiss 
law x(T) = C/(T + V2T£) with T£ rs 35 K and an 
effective paramagnetic moment of 2.5 /Xb- In contrast, 
at low T, x(T) increases with La dilution, reaching the 
single-impurity limit for x > 0.95. Similar systematics 
are also observed for the ab-plane component, Xab(T). 

We performed a CEF fit to the ab-plane and c-axis 
components of the susceptibility for the dilute Ce regime 
(see solid line in the inset of Fig. 2). The ground state is 
found to bealf doublet (— 0.39| ± 5/2) + 0.92| q= 3/2), 
including a low-temperature Weiss temperature of 2 K) 

(2) 

with a fairly large gap to the excited Ty ^ (0.92 1 ± 5/2) + 
0.39| T 3/2), Ai = 148.0 K) and r 6 (| ± 1/2), A 2 = 196.5 
K) doublets. This scheme also well fits the high T sus- 
ceptibility above 80 K for all x for both the ab-plane and 
c-axis components, which suggests that the CEF scheme 
is essentially independent of x. In order to confirm the 
CEF scheme, we measured the specific heat up to 200 
K for x = 0, 0.25, 0.50, 0.75, 1.00. The magnetic con- 
tribution to the specific heat C m is estimated as Cp of 
Cei-^La^CoIns minus Cp of LaCoIns divided by the Ce 
concentration. For all x, C m shows a broad peak around 
100 K, consistent with the Schottky type fit based on the 
above CEF scheme (see Fig. 3 right inset). These data 
also imply that the CEF scheme is roughly independent 
of x. There are differing CEF fits for CeCoIns in the lit- 
erature. Our results agree best with the scheme of Rcf. 
JTo| | and less well with those in Refs. |l3||. We argue that 
these discrepancies arise from the way intersite correla- 
tions are taken into account. 
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FIG. 3. C m /T times T* s vs. T/T* a . The solid line repre- 
sents the — InT fit. Left inset: T dependence of C m for x — 
0.98 and 0.99. The solid curve is the fit to the S = 1/2 Kondo 
impurity limit with Tk = 1.7 K. Right inset: T dependence 
of C m . The solid curve is the fit based on our CEF scheme. 



Given the essential same CEF scheme for all x, we in- 
terpret the change in x(T) with La dilution at low T 
as an intersite coupling effect. To quantify this contri- 
bution, we subtract the single impurity CEF x(T) from 
x(T) for each x to get the "excess susceptibility" Ax(T), 
which shows a strong decrease of % at a concentration de- 
pendent characteristic temperature T* (e.g., around 50 
K for CeCoIn 5 ). Interestingly, we found a scaling law 
for both ah and c components of A%(T) vs. T/T*(x) as 
shown in Fig. 2. Basically the same systematic change 
of T* with x is found for both components of %, T* ab and 
T* c , as seen in Fig. 4. 

For Kondo systems C m /T is given by 1/T* S (T* s 
is a characteristic energy) times a scaling function of 
T/T* s [|. Actually, we have found two different scal- 
ing functions in Cei_ K La x CoIn5. In the coherent re- 
gion (a: < 0.5), for the temperature range between 
0.04T* S and 0.6T* S , C m /T fits well to the expression 
-5.4/T* s ln(T/T* s ) (J/mol-Ce K 2 ) as displayed by the 
solid line in Fig. 3. Considering 5.4 J/mol-Ce K is close 
to R In 2, this expression suggests that T* s gives the tem- 
perature scale to recover the ground doublet entropy. In 
the incoherent regime (x > 0.5), all the curves scale to 
another function as shown in Fig. 3. Furthermore, the 
C m /T data for x — 0.98 and 0.99 agrees well with the ex- 
act results for the S — 1/2 single- ion Kondo model with 
the single parameter Tk = 1.7 K, which defines the T* s 
at the single impurity limit [|l4| (see left inset of Fig. 3). 

Generally, Tk is related to the high-temperature T% 
and the CEF splittings via k B T K = (fc B T£) 3 /Ai A 2 
From the susceptibility we have T^ as 35 K, Ai = 148 
K, and A 2 = 196.5 K, such that the expression yields 
T K ~ 1.5 K, consistent with the T* s and the low T Weiss 
temperature in the single impurity limit. 
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FIG. 4. Energy scale diagram of Cei-zLa^CoIns. Tsc is 
the superconducting transition temperature determined by 
the onset of the jump observed in Cp(T). Tii near and Tfl are 
the temperatures below which the resistivity starts to show T 
linear and T 2 dependence. T* ab , T* c and T* s are the charac- 
teristic temperatures of the ab-plane and c-axis susceptibili- 
ties, and the specific heat, respectively. 

With the scaling procedures above, we obtain a sys- 
tematic La concentration dependence of T* s as shown 
in Fig. 4. There are three points to note. First, the 
characteristic temperatures T* s , T* ab and T* c are essen- 
tially identical, indicating that there is only one energy 
scale besides Tk. Second, in the x — > 1 limit this energy 
scale agrees with Tk, showing that the scale T* originates 
from the single ion Tk- Third, in the dense Kondo limit 
CeCoIns, T* coincides with the T co h of the resistivity. 

Because Ig is constant throughout the alloy series, and 
the high T \ an d Cp data point to no changes in the CEF 
parameters, the single ion Tk should also be constant. 
The systematic increase in T* must then arise from the 
intersite coupling. Hence, this system has three well- 
separated energy scales: Tk ~ 1 — 2 K for the single ion 
Kondo screening, T* w 40—50 K for Kondo lattice effects 
and the CEF excitations A w 200 K. The fact that these 
energy scales are very different is a fortuitous situation 
that allows the present scaling analysis. 

The different scaling of C m /T in the coherent and in- 
coherent regimes suggests that the nature of the ground 
state changes near x = 0.5. A qualitative change was 
also found in the low-T dependence of the resistivity. 
The right inset of Fig. 1 shows the log-log plot of the 
inelastic part of the in-plane resistivity vs. T/T co h in the 
coherent regime. All curves for x < 0.25 have a T-lincar 
dependence below a characteristic temperature Tii ncar . 
Note that all the T-linear curves have the same coeffi- 
cient A/T co h, which reveals that the inelastic scattering 
rate scales with the l/T co h. Especially for CeCoIns, it 
scales with 1/T* because T co h = T* at x — ► limit. 

In contrast, the left inset of Fig. 1 displays the T 2 
dependence for x = 0.75 and 0.955 in the incoherent 
regime, which becomes evident below the Fermi liquid 
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temperature Tfl — 100 mK. In the single impurity limit, 
according to Nozieres and Yamada Jl(J the T 2 coefficient 
of the conductivity should be proportional to T^ 2 . Using 
their formula, 7k is estimated to be 1.1 K for x = 0.955. 
This is consistent with the estimates above. 

Non-Fermi-liquid (NFL) behavior is manifested in the 
T-linear resistivity and - hx(T/T*) dependence of C m /T 
(Figs. 1 and 3). NFL behavior has been observed in 
numerous other Ce, Yb and U alloys and compounds 
jl7j and is frequently attributed to a quantum critical 
point (QCP) due to the proximity of an AF instability. 
Especially near a QCP of 2D antifcrromagnetism, criti- 
cal fluctuations are expected to lead to p cx T/T* and 
C m /T cx - \n(T/T*) as experimentally observed. 

In this case, T* gives the energy scale of the spin fluctu- 
ations jl|-|ij. 

In addition, several experiments actually point to the 
2D AF critical fluctuations. De Haas-van Alphen mea- 
surements have shown that Cei_ x La x Coins has a strong 
2D anisot ropy o f the Fermi surfaces reflecting its layered 
structure JK J2G] . The resistivity and specific heat under 
pressure |lC| ll"|] and NMR 1/T\ measurements [^l| sug- 
gest that CeCoIns is close to the QCP of quasi-2D anti- 
ferromagnetism. Optical measurements of CeCoIns have 
revealed the formation of an additional peak below 50 K 
within the hybridization gap 22 1, which is attributable 
to the quasi-particles interacting with AF fluctuations. 
Notably, their estimate of the fluctuation energy yields 
the same energy scale (8 meV) as T* for CeCoIns (45 K). 

Moreover, the decrease of the excess susceptibility A% 
below T* is almost isotropic as shown in Fig. 2, and sug- 
gests the development of RVB-like singlets (short-range 
AF correlations of / moments). The A\(T) is well- 
described by the expression Ax(0)[l — exp(— 0.05T*/T)] 
with Axab(O) « -0.038 emu/mol-Ce and Ax c (0) ~ 
—0.052 emu/mol-Ce. The reduced activation energy 
A = 0.05T* suggests slightly dispersive AF correlations. 
Interestingly, the A for CeCoIns gives the same energy 
scale as the possible psuedo gap temperature (w 3 K) 
where the resistivity decrease was observed pd| . 

These results suggest that spin liquid formation with 
the gap A starts through the screening process, involv- 
ing the intersite AF correlations with the energy scale T* . 
Furthermore, the qualitative change in the T dependence 
of both C m /T and p m around x = 0.5 indicates that a Ce 
concentration larger than 50 % is the criterium not only 
for coherent transport but also for the intersite coupling 
leading to a short-range collective mode. This implies 
that the correlation length is of the order of the lattice 
constant a and is rather short compared to what is usu- 
ally expected for the RKKY interaction. 

In summary, we have studied the effects of La dilu- 
tion on CeCoIns. The scaling laws found in the sus- 
ceptibility and specific heat reveal a systematic evolu- 
tion of the characteristic intersite coupling energy T* . 
Below this temperature, we found a Fermi liquid state 
in the incoherent impurity regime, while NFL behavior, 
which scales with T* has been observed in the coherent 



dense Kondo regime. The latter suggests that the Kondo 
screening in the lattice is no longer local but rather an 
intersite effect involving AF spin fluctuations. Given 
that the ratio of Tsc to the spin fluctuation energy T* 
for CeCoIns (ps 5 %) is of the same order as those for 
cuprates and other heavy fermion superconductors [fTs] , 
it is suggestive to associate T* with the energy scale lead- 
ing to the relatively high Tsc of CeCoIn 5 . 
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